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Anomalous Hall effect driven by dipolar spin waves in uniform ferromagnets 
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A new type of anomalous Hall effect is shown to arise from the interaction of conduction electrons 
with dipolar spin waves in ferromagnets. This effect exists even in homogeneous ferromagnets with¬ 
out relativistic spin-orbit coupling. The leading contribution to the Hall conductivity is proportional 
to the chiral spin correlation of dynamical spin textures and is physically understood in terms of 
the skew scattering by dipolar magnons. 


The anomalous Hall effect (AHE) has been known for 
more than a century and provided various intriguing ar¬ 
eas of researches both theoretically and experimentally 
[T] . The AHE are due to both the presence of the sponta¬ 
neous magnetization and spin-orbit interaction, and thus 
no external magnetic field is needed as opposed to the 
ordinary Hall effect. The contributions to the AHE are 
conventionally understood by the intrinsic and extrinsic 
mechanisms; the former is associated with the anomalous 
velocity deriving from non-trivial band curvature (Berry 
curvature), and the latter comes from spin-orbit scat¬ 
tering from impurities [2]. Distinctly from these mecha¬ 
nisms, another type of contribution to the AHE, known 
as the topological Hall effect (THE), was later proposed 
[311]. The THE arises from chiral magnetic textures such 
as Skyrmion configuration, which generate Berry phase 
for conduction electrons leading to Hall effect [5|- This 
Berry phase interpretation is appropriate in the regime of 
strong coupling between conduction electrons and mag¬ 
netic textures. On the other hand, the Hall conductivity 
in the weak coupling regime was shown to be proportional 
to the scalar spin chirality, given by ■ {Sr2 x Srs), 
which measures the correlation between the spin orien¬ 
tations and real-space positions [B|- The leading order 
contribution reads 

CC ^ ' {^r2 X ^ z ^ 5 (1) 

ri,r 2 ,r 3 

where a = ri — r2, b = r2 — r^, c = — ri, T (a, b, c) 

is a symmetric function of a = |a|, b = |b| and c = |c| 
related to electrons, and the summation is taken over all 
locations of localized magnetic moments ■ In order 
to have a non-vanishing contribution in this setting, the 
scalar spin chirality must include a component that is 
odd under parity (i.e. n = {xi,yi,Zi) (ix*, =F2/i, 
which, in translationally and rotationally symmetric sys¬ 
tems, requires ■ {Sr2 x S',.3)|odd oc (a x b)^. 

Recently, Hall effects not only of electrons but also of 
various bosonic excitations, involving photons, magnons 


and phonons, have been studied |7]. Amongst them is 
the Hall effect of magnetostatic forward volume waves, 
or dipolar magnons, which have Berry curvature origi¬ 
nating from the magnetic dipole-dipole interaction |5]. 
Since conduction electrons interact with magnons via the 
exchange interaction, it seems natural to expect that the 
electrons may also acquire some chirality in their orbital 
motion and exhibit a Hall effect. 

In this Letter, we approach this problem by a pertur¬ 
bative expansion in the exchange coupling between the 
dynamical magnetization and conduction electrons. We 
derive a spin-chirality formula for the Hall conductiv¬ 
ity, which can be regarded as a generalization of Eq. 0 
to spin texture fluctuations. It is found that the dipole- 
dipole interaction induces a parity-odd component of chi¬ 
ral spin correlation, 

{Sri ■ iSr2 X Sr^)) |odd = (a X 6)^ (a • 6)j F{a)F{b) 

-bcyclic perms in {a,b,c), (2) 

where (■•■) denotes quantum/thermal averaging, the 
subscript t indicates projection onto x-y plane (a • b)t = 
Qxbx + ciyby, and F{a) is a function related to magnons. 
Because of the factor (a • b)t in Eq. 0 , there is no net 
spin chirality in the sense that the spatial average of 
{Sri ■ i^r2 X 5^3)) (a X 6)2 is zero. Nevertheless, the 
electron contribution, which corresponds to (a x b)zJ^ in 
Eq. ([^, provides factors of (a • b)t in such a way to ren¬ 
der the Hall conductivity finite, as we shall demonstrate 
below. 

Let us first reformulate the theory of dipolar magnons 
in a form suitable for studying their interaction with elec¬ 
trons. We consider a normally magnetized film of ho¬ 
mogeneous ferromagnet. The magnetization rfi (r) satis¬ 
fies (rh (t'))q = (0,0, Mg) in the ground state, where Mg 
is the saturation magnetization. Let us introduce the 
Holstein-Primakoff field Oq by 



q 
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where m± = rhx ± irhy, 7 < 0 is the effective gyromag- 
netic ratio of electron, and V is the volume of the sample. 
In the presence of both exchange and dipole-dipole inter¬ 
actions, the magnon Hamiltonian can be written as [5] 

= ( 4 ) 

Q 

in terms of the diagonalizing operator, bq = 

A~e^‘^'’a^_q, and the eigenfrequency 



where q = q{sm9qCOS<j)q,sin9qSm(j)q,cos9q). The Bo- 
goliubov coefficients are given by 

_ gjq ± (cup -f 

2^/uiq {ujo+ujmM'^) 

If we neglect magnetocrystalline anisotropy, ujq = —yiJex 
and lom = —iMg, where i7ex is an external magnetic 
field, and A is the square of the exchange length, whose 
typical value is ~ 10“^® [m^]. Note that Og is not an eigen 
operator unless the dipole-dipole interaction is absent, 
Aq =0. To treat interactions between magnons and 
electrons, we shall apply quantum field theory techniques 
[BE] and define thermal Green’s functions 

Vab (1,2) = (T^ [rhA (Ti,ri)mB (r2,r2)]) , (7) 

where T,- indicates imaginary-time ordering, A, B = zt, 
and mAir) = A consequence of the dipole- 

dipole interaction is the spin-flip propagator 

2?±± = - ^ E , (8) 

l.q 

where vi is the Matsubara frequency and Vq (ivi) = 
—2tujjql{v‘j -I- fi^ojq). We neglect the self-energy for sim¬ 
plicity. The phase factor in Eq. (|^ is a manifes¬ 

tation of the coupling between the spin and orbital mo¬ 
tion of magnons. In particular, the spin-flip propagator 
Q contains a parity-odd element, which will ultimately 
lead to electron Hall effect. The spin-preserving part 
I?±zp also receives a contribution from the dipole-dipole 
interaction, but it is parity even and does not essentially 
affect AHE. 

For the electron part of Hamiltonian, we consider a 
free electron gas with the s-d exchange coupling to rh; 

_^ 

Hs = 'y ) ^ 

k 

■ (4'^Cfe+g) " . (10) 

Here, = (cj[,.^,c^) is the electron creation operator, /i 
is the chemical potential, J the s-d exchange energy, N 



FIG. 1. (color online) Chiral spin correlation x|odd. (a) As a 
function of ra = {x, y, 0 ) for ri = ( 0 , 200 , 0 ), ra = ( 0 , — 200 , 0 ) 
and ri = r 2 = T 3 . The unit of length is A. (b) As a function 
of a = |ri — ral at {x,y) = (100,0), (100,400) and (400,400) 
plotted by red, blue and green lines, respectively. As seen, 
xiodd is odd in x while oscillating in x,y and a, consistent 
with the vanishing spatial average of (a x ^)zX- 

the total number of localized magnetic moments in the 
sample and cr the triplet of the Pauli matrices. Our goal 
is to compute the Hall conductivity a^y based on the 
Hamiltonian H = Hs-\- Hgd and the Kubo formula 

cTxy (w) = — f dr (T^ [jx (t) jy]) , (11) 

lUJ J 

where j = is the charge-current density 

operator, and the analytic continuation iuj\ uj -\- iQ 
is understood. The electron propagator is given by 
Qku (f^n) — ^ka (fs^n)} with ^ka — Ck “t” 

JSa, where S = —MsV/h'jN is the magnitude of lo¬ 
calized spins, a = ±1 the spin projection, and E is the 
self-energy. We take ^/.^(e-l-iO) = —frfccr(e), where T is 
the damping parameter, neglecting its real part. 

We assume J is small and calculate <Jxy perturbatively 
with respect to J. As in Ref. [6], a non-trivial contribu¬ 
tion shows up at the third order. The spin dependence of 
Qkcr turned out to be unimportant (not essential to AHE, 
just giving higher-order corrections), which will thus be 
neglected in the following: = Qk- Using the identity 

ir {uiOjCTk) = 2ieijk, we obtain 

(Tr [jx (t) jyj) = J 

X 2i (T, [rh (1) • {m (2) x m (3)}]) £xy (1, 2, 3) , (12) 

where £xy (1,2,3) is the electron part consisting of t/’s, 
and 1 = (ri,Ti), J dl = f dri J dri, etc. This expression 
can be regarded as a generalization of the spin-chirality 
formula of Hall conductivity obtained in Ref. |H] to the 
situations where the spin texture is allowed to fluctuate. 

In order to have non-vanishing axy, it is necessary that 
the time-ordered chiral spin correlation 

X = (T^ [rh(l) • {rh(2) X rh(3)}]) (13) 

has parity-odd components since £xy (1,2,3) is parity 
odd. As anticipated, the spin-flip magnon propagators 
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FIG. 2. The diagrams contributing to a:^y 


at J^. In terms of 
£xy (1, 2 , 3), there are only two configurations; those in the top 
row and those in the bottom. The three variations within each 
row occnr due to the cyclic permutation of magnon arguments 
in Eq. (14l. Electron-magnon vertices with one (two) wavy 
line(s) represent x, y components {z component) in Had- 


of Eq. ([^ serve this purpose. In fact, one can derive 
from Wick’s theorem the parity-odd part of x 8 -s 


X lodd = (1, 2) V++ (1, 3) (14) 

— I?++ (1, 2) V _(1, 3) -I- (cyclic perms) . 


Here we expanded to the lowest order in the magne¬ 
tization fluctuation as m-z = Mg — TO_m+/(2Ms). More 
explicitly, this takes the form 


^ lodd 


{a X b)z 


{a ■ b)tFaFb + {b ■ c)tFbFc 


+ (c • a)tFcFa , 


(15) 


where Fa = (oa, ± tOy) ^'D±±{a-,Ta) depends only on 
at = \ + a^^az and Ta = tx — T 2 (i.e. Fa is a scalar 


with respect to rotation around z-axis). Equation (15) is 
an elaboration of Eq. (^ 2 ^. Spatial variation of y I jj is vi- 


odd 


sualized in FIG. We used A = 0.4 x 10“^® [m^], liwo = 
hujM = 10“^[eV],fcBr = 10“^ [eV] as a representative 
set of parameters. The q-integral in Eq. (|^ was evalu¬ 
ated by introducing a cutoff gm = 1 [A ]. Remarkably, 
there is a definite chiral spin correlation with odd parity 
once the relative locations (a and b) are fixed. If aver¬ 
aged over the locations, however, it vanishes because of 
the factor (a ■ b)t, as mentioned in the introduction. 

The Hall conductivity a^y is obtained by combining 
the magnon part y with the electron part £xy (1,2,3). 
The resulting diagrams are shown in FIG. where the 
solid and wavy lines represent Q and V, respectively. Ev¬ 
ery term of £xy is a product of a factor (a x b)z (or its 
equivalent) and functions of a, b and c. Since c = |a -|- 6 | 
contains a factor (a • b)t, the sign change due to (a ■ b)t in 
Eq. ([^ upon integration is compensated, resulting in a 
non-vanishing cr^y. This is to be compared with the case 
of Ref. [6] , in which the electron part that is independent 
of (a • b)t contributed to cr^y. 

To evaluate axy explicitly, it is convenient to work 
with the Fourier (wave-vector and Matsubara frequency) 
space. After performing analytic continuation zwa —>■ 
hoj + iO and taking the limit w —>■ 0 , one finds a number 
of terms containing different combinations of retarded 
and advanced Green’s functions. As we are computing 
a transport quantity, these terms are naturally classified 
according to the powers of the electron lifetime, or de¬ 
cay rate T^ (e). In the clean limit, the dominant term 
is that of the lowest power in T [ 1 ], originating from the 
two pairs of retarded and advanced electron propagators 
at both ends of the diagrams with identical momentum 
arguments. This is similar to the so-called extrinsic con¬ 
tribution and reads 


mu — 


2h^e'^J^SV^ 


J de J dei j *2 ^ A+sin [2 (())q - ())q/)]/'(e) 

k-q.q' 


X (g - e) + / (g)} {n (£2 -e)+f (£2)} (£)| ImH^ (£1 - e) ImD^, (€2 - e) 

{(k + g) X (k + q')}JmG^(€) (£i)|^ - 2 (k x g')^ | G^ (£) | ^ ImG^+g (£1) 


(16) 


where n(x) and f(x) are Bose and Fermi distribution 
functions, and G^ and are the retarded counterparts 
of Q and (D, respectively. In the square bracket in the 
last line of Eq. (16 1 , the first term comes from the di¬ 
agrams in the first column of FIG. while the second 
term from the sum of the rest. To the lowest order in T, 
one can replace [G^ (e)] by nS (e — ^k) /Efc (e) and other 
propagators by their unperturbed form. The products of 
the distribution functions effectively restrict the energy 


variables £, £1 and £2 to [—/csT, fc^T], much smaller than 
and the delta functions force k, k + g and k + g' to 
be near the Fermi surface and T’s can be regarded as a 
constant Tp- = rfcj,,(ei;’). Integrating over energies and 
(j)k, we arrive at 


POO POO 

o'xy = crojl^ / dkz / dktf'(^k)kt (17) 

J — 00 Jo 

X{lx{k)-h{k)}{h{k)-2l3{k)} , 
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FIG. 3. (color online) (a) (7i - l 3 ){Ii - 2/3) in Eq. 
for (ktykz) = fcF(sin0,cos0) as a function of cos9. Note 
that ktdO = —kFd{cos6). (b) Hall conductivity evaluated 
by Monte Carlo integration (red dots) with a fit \axy\/oQ = 
6.83 X 10“® X y 2 . 29 g-o.ooi/T Plain power-law fit¬ 

ting gives ~ j'2'48 qijjg same parameter set as in FIG. has 
been used along with /r = 0.6 [eV]. 


where CTo = Se'^kp/{‘^irheF^^), k now denotes doublet 
{kt^ fcj), kt = ^Jk‘^-\- ky and the functions (fc) are given 

by 

^ ^ V t (Pq ujm 
" “ ^2mN ( 2^)3 ivq 

X {n (hWq) + f (a^k + kujq)} 6 {^k+q - ■Cfc - akuq). 


{k ■ q)t 


qth 


Qt 


(18) 


We note a^y is negative for JS > 0. This is seen by ap¬ 
proximating —/' by delta function; then the integral is es¬ 
sentially the area under the curve in FIG.[^(a). Temper¬ 
ature dependence has been studied numerically, yielding 
a roughly power-law dependence Oxy oc with a small, 
but appreciable, deviation [FIG. (b)]. The estimated 
value of p falls in the range 2.2 < p < 2.5 depending on 
the choice of fitting function. 

The mechanism of the present AHE can be identified 
as the skew scattering by the chiral spin correlations of 
dipolar magnons. In fact, by relabelling the wave-vectors 
in FIG. the on-shell scattering amplitude of two elec¬ 
trons with momenta p, p' at the ends of the diagrams 
going through y [odd (with an intermediate electron line) 
can be written as proportional to (p x p') 2 . 

To conclude, we have shown that there is a nonzero 
contribution to the AHE originating from the interac¬ 
tion with the dipolar magnons in homogeneous ferromag- 
nets. The leading order term occurs at the third order in 
the s-d exchange interaction and involves a factor of the 
chiral spin correlation felt by traversing electrons. The 
result serves as a generalization of the Hall effect from 
static spin chirality as well as giving a concrete example 
in which anomalous Hall current requires neither rela¬ 
tivistic spin-orbit interaction nor fixed magnetic texture. 
In a broader sense, the dipole-dipole interaction acts as a 
kind of spin-orbit interaction in that it couples the spin 
and orbital motion of magnons. Electrons then inherit 
the chiral orbital motion through skew scatterings by the 
quantum/thermal dipolar magnons. 


Unfortunately the predicted Hall conductivity is small 
because of the dominance of exchange energy over the 
dipole-dipole term, i.e. \q^ ^ 1 for typical < 7 , which is 
y/kBT/{hu;M^)- To identify the proposed AHE experi¬ 
mentally, therefore, material should be chosen carefully. 
Eirst of all, conventional AHE due to spin-orbit interac¬ 
tion must be suppressed in conductors; for example, in 
some ferromagnetic metals, anomalous Hall coefficients 
are known to be compensated by blending elements with 
opposite signs of AHE. Since further numerical studies 
(not shown) indicate that jcrxyl decreases with as long 
as ep ^ kpT, materials with small Eermi energy will 
be favorable, such as a semiconductor magnet or a hy¬ 
brid film comprising a ferromagnetic insulator layer and 
a low-carrier conductor layer. In many materials, AHE 
accompanies spin Hall effect which is often used in spin- 
tronics, and the present mechanism may open the door 
to spintronics free from spin-orbit interaction. 
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